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ABSTRACT 

We present a mean field model based on the approach taken bv lRempel <|2005|> in order to investigate 
the influence of stochastic fluctuations in the Reynolds stresses on meridional flow and differential 
rotation. The stochastic fluctuations found in the meridional flow pattern directly resemble the 
stochastic fluctuations of the Reynolds stresses, while the stochastic fluctuations in the differential 
rotation are smaller by almost two orders of magnitude. It is further found that the correlation length 
and time scale of the stochastic fluctuations have only a weak influence on meridional flow, but a 
significant influence on the magnitude of variations in the differential rotation. We analyze the energy 
fluxes within the model to estimate time scales for the replenishment of differential rotation and 
meridional flow. We find that the time scale for the replenishment of differential rotation (~ 10 years) 
is nearly four orders of magnitude longer than the time scale for the replenishment of meridional flow, 
which explains the differences in the response to stochastic fluctuations of the Reynolds stress found 
for both flow fields. 

Subject headings: Sun: interior — rotation 



1. INTRODUCTION 

Recently iRempell l|2005f) presented a mean field 
model for solar differential rotation and meridional 
flow, assuming a parameterization of all convective 
scale processes, most importantly the turbulent an- 
gular momentum transport (A-effect) introduced by 
I Kitchatinov fc Riidigerl l) 19931) . The model presented by 
IRempell 1)20051) and also ot her mean field models for dif- 
ferential rotation (see e.g. IKitchatinov fe Rudigeilll995t 
iRiidiger et al.lll998HKiiker fc Stixll2001[) use a time inde- 
pendent parameterization of the turbulent angular mo- 
mentum transport leading to stationary solutions de- 
scribing the mean flows. On the other hand, 3D nu- 
merical simulations show a significant temporal varia- 
tion of meridional flow and differential rotation due to 
the fluctuations in the convective motions leading to the 
Reyn olds stresses ( Micsch et alJ 120001 iBrun fc Toomrd 
2002). 

Observations of differential rotation are usually based 
on 3 to 4 month sets of data and therefore can only tell 
us somethin g about variations o n time scales longer than 
this period l)Schou et alJ 120021) . The la error-intervals 
of inversions within the convection zone are typically 
around a few nHz or less than 1% of the rotation rate. 
The differential rotation also shows a systematic cycle 
variation, the torsional oscillations, which have an am- 
plitude of around 1% of the rotation rate. By contrast 
the variability found in observations of the meridional 
surface flows is much larger compare d to the average am- 
plitud e. Over time scales of years iZhao fc Kosovichevl 
( 2004) found variations of the meridional surface flow of 
the order of 10ms -1 compared to a mean flow amplitude 
of 20ms -1 . While these variations are most probably also 
related to the solar cycle, much less is known about short 
term variations, which could be related to the turbulent 
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origin of the flow itself. Comparison of the long term 
variations found in differential rotation and meridional 
flow indicates that the relative amplitude of the merid- 
ional flow variability is almost a factor of 100 larger than 
the variability of the differential rotation. 

In this paper we present a mean field model for differen- 
tial rotatio n and me ridional flow following the approach 
described in lRemneH l|2005j) to investigate the influence of 
stochastic fluctuations in the Reynolds stresses on merid- 
ional flow and differential rotation. The main intention 
of this paper is to investigate how much variability can 
be expected in the meridional flow pattern, given the 
strong constraints set by helioseismology on the variabil- 
ity of the 3 to 4 month mean of the differential rotation. 
The possible influence of the solar cycle on both flow 
patterns will addressed later in a separate paper. 

2. MODEL 

For thi s investi g ation we use the model described in 
detail in IRempell l)2005l) . We add random noise with 
defined correlation length and time scales to the param- 
eterization of the turbulent Reynolds stress responsible 
for driving the diff erential rotation and meridional flow 
[see IRempell l)2005|) . equations (31) and (32)]: 

A r0 = A 0r = +L(r, 9) cos(6» + A(r, $)) 

(l + cCr(r,e)/a r ) (1) 

A e0 = A 0e = -L(r, 9) sin{9 + A(r, 9)) 

(l + c( e (r,e)/ae) (2) 

L(r, 9) denotes the mean amplitude of the turbulent an- 
gular momentum flux, whereas A(r, 9) describes the mean 
inclination of the flux vector with respect to the axis of 
rotation. L(r, 9) and A(r, 9) need to be antisymmetric 
across the equator to fulfill the sym metry constrain ts of 
the A-effect. For further details, see IRempell l)2005j) . 

£ r and £g denote random functions, which will be de- 
scribed below in more detail. The parameter c deter- 
mines the amplitude of the random noise in units of the 
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Fig. 1. — Variance of meridional flow (left panels) and differential rotation (right panels) at the top of the domain r = O.985i?0. The top 
panels a) and b) show the total variance of the time series for the cases 1 to 4 (variation of correlation time scale) . The solid line indicates 
case. 1 (r = 0.25^0 1 ), the dotted line case 2 (r = 1Q 1 ), the dashed line case 3 (t = 4f2 Q ), and dashed-dotted line 4 (t = 16f2 Q ). 
The middle panels c) and d) show cases 5, 3, and 6 (variation of correlation length scale). The solid line indicates case 5 (A = 0.05), the 
dotted line case 3 (A = 0.1), and the dashed line case 6 (A = 0.2). The bottom panels e) and f) are similar to the top panels, however the 
variances are computed from a four month mean. 



standard deviations a r and ag. By perturbing both com- 
ponents of the turbulent angular momentum flux with 
uncorrelated random functions we allow for a change of 
amplitude and direction of the angular momentum flux. 

In order to generate a random function with a defined 
correlation length scale, we construct a 2D random field 
in the r-0-plane by superposing Gaussian functions with 
a fixed width that determines the length scale. The po- 
sition in the r-#-plane of individual Gaussian as well as 
the amplitudes are random. We introduce a correlation 
time scale into the problem by performing a temporal 
average over the 2D random fields. In principle this av- 
eraging could be done by creating a long time series of a 



2D random field and using then a running mean. How- 
ever, this approach would require a fairly large amount 
of memory. Instead we produce a 2D random field g each 
time step and solve the following equation to introduce 
a correlations time scale: 

Ot T c T c 

where £ denotes the random field with a correlation time 
scale r c . The main disadvantage of this approach is that 
the temporal average is always dominated by the most 
recent realizations of g. This effect can be reduced by 
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TABLE 1 

Parameters for the 
stochastic 
perturbations 



case t c [Qq 1 ] 


A 


1 0.25 


0.1 


2 1.0 


0.1 


3 4.0 


0.1 


4 16.0 


0.1 


5 4.0 


0.05 


6 4.0 


0.2 


Note. 


Models 



1 to 4 cover a range 
in the correlation time 
scale by a factor of 64, 
whereas models 5, 3, 
and 6 cover correlation 
lengths scales varying by 
a factor of 4. Qo denotes 
the rotation rate of the 
solar interior. 



applying this method recursively as follows: 
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T c /n 
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T c /n 


T c /n 
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- 9 , 


dt 


T c /n 


T c /n 



(4) 



where £i is the random field used for the simulation. The 
functions C2 to ( n are auxiliary functions required for the 
averaging process. The results presented in this paper are 
obtained using a value of n = 5 for the temporal averag- 
ing. The free parameters in this approach for generating 
a 2D random field are the correlation length scale A and 
the correlation time scale r c . A measures the half width 
of the Gaussians at half maximum in units of the domain 
size in the radial and latitudinal direction (therefore ra- 
dial and latitudinal width are not the same). For the 
results presented here we use an amplitude [defined by 
the parameter c in equations and of 0.33, mean- 
ing the stochastic fluctuation imposed on the Reynolds 
stress has a la fluctuation corresponding to 33% of the 
mean. 

In order to allow for a detailed statistical analysis of 
meridional flow and differential rotation, we compute 
during the simulation averages of all variables and the 
square of all variables over the output sampling period, 
which allows computing the mean flows and the stan- 
dard deviations including all simulation time steps and 
therefore avoiding the influence of the output sampling. 

3. RESULTS 

For evaluating the effect of random noise in the tur- 
bulent angular momentum transpo r t we u se the model 
described in the case 1 of iRempell l)2005l) as reference. 
Table Q summarizes the parameters we use for the gen- 
eration of the random noise fields used to randomize the 
Reynolds stress according to equations JTJ and All 
other model parameters are given in lRempell (2005). The 



cases 1 to 4 cover a range in the correlation time scale 
by a factor of 64 (from roughly one day to two months), 
while the series of cases 5,3, and 6 encompasses a varia- 
tion in the correlation length scale of a factor of 4 (from 
a radial half width of about 10 Mm to 40 Mm). 

In the following discussion we will focus on the variance 
of the meridional flow a v and the variance of the differ- 
ential rotation an- Since helioseismic measurements of 
differential rotation typically involve 4 month averages, 
we also compute the variance of 4 month averages of our 
model output. 

FigureQshows a v (left column) and an (right column) 
computed at the top of the domain (r = O.985i?0) for the 
models 1 to 6 listed in Tabled The top panels show the 
total variability for the cases 1 to 4 (change of correlation 
time scale), while the middle panels show the cases 5, 3, 
and 6 (variation of correlation length scale). The bottom 
panels show the same properties as the top panels, but 
computed from a four month mean. 

The most striking feature visible in panels a) to d) is 
that while the fluctuations of the meridional flow show 
no systematic variation with correlation time and length, 
the fluctuations of differential rotation show a systematic 
increase with both. An exception is case 2 in panel a) and 
b), since r c = Q^ 1 is very close to the time scale of iner- 
tial oscillations and therefore a v is enhanced because of 
a resonance. The amplitude of a v is about 20% to 25% 
of the mean meridional flow (in our model the merid- 
ional flow reaches a peak value of around 13.5 ms _1 at 
45° latitude), which is close to the noise level added to 
the Reynolds stress (33%). On the other hand, except 
for high latitudes the value of an is less than 1% of the 
reference rotation rate in all cases. 

Panels e) and f) show the similar quantities as panels 
a) and b) , however computed from a four month mean of 
the time series, similar to the averaging interval typically 
used in GONG/MDI inversions of differential rotation 
profiles. It is not surprising that the amplitude of the 
fluctuations decreases, however the effect is stronger in 
case of the meridional flow, e.g. the averaging decreases 
a v for case 3 (dashed line) by a factor of 2, while an is 
nearly unaffected. 

So far we focused on the variability of meridional flow 
and differential rotation at the top of the domain. Fig- 
ure [3 shows the depth dependence for the meridional 
flow at 45° and the for differential rotation at the se- 
lected latitudes 0°, 15°, 30°, 45°, 60°, and 90°. In both 
cases solid lines indicate the mean and the grey shade 
the la interval. Shown is case 3 with a correlation time 
scale of 4J7 ~ 1 . Both a v and an show a maximum at the 
top of the domain and a nearly monotonic decrease to- 
wards the base of the convection zone. We emphasize 
that this result was obtained by using a constant corre- 
lation time and length scale throughout the convection 
zone. In a more realistic description it can be expected 
that both correlation time and length scale decrease to- 
wards the solar surface, which according to the trends 
found in Figure would lead to a decrease of an close 
to the surface, while a v remains more or less unchanged 
(except for the enhancement of the signal close to the 
resonance). In our investigation we also use a constant 
amplitude of the stochastic fluctuations of the Reynolds 
stress throughout the convection zone. A depth depen- 
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Fig. 2. — Depth variation of the standard variations of meridional flow and differential rotation for case 3. Panel a) shows the meridional 
flow profile at 45° latitude and the la intervals as grey shade. Panel b) shows the differential rotation at 0°, 15°, 30°, 45°, 60°, and 90° 
latitude with lcr intervals as grey shade. 



dence of this amplitude would influence both, a v and ctq 
in a similar way. 

We computed all models shown in this paper with an 
amplitude of c = 0.33 for the random noise in the A- 
effect. Note that a v and ctq scale linearly with c, meaning 
that the results given here can be rescaled to apply to a 
different amplitude of the random noise. 

We indicated already above that the variance of the 
meridional flow and differential rotation is enhanced 
around r c ~ fl^ 1 because of the excitation of inertia 
oscillations. Figure shows the variances of meridional 
flow and differential rotation at the top of the domain 
(r = 0.985 R Q ) for the latitudes 22.5° (solid), 45° (dot- 
ted), and 67.5° (dashed) as function of r c . Panel a) shows 
the total variance of the meridional flow, clearly indi- 
cating a resonance around r c ~ f — 2 . Since our 
time scale r c is defined as an e-folding time scale for the 
random perturbations, we expect this resonance around 
t c w Ti/2 = n/uJi, with the frequency of inertia oscilla- 
tions given by uji — 2f2o- This leads to a locations of the 
resonance around 

< 5 > 

Figure shows a weak tendency of increasing r r with 
decreasing latitude, which indicates that the geometry of 
the problem and the stratification leads to some confine- 
ment of motions in latitude (for motions in latitude only 
we would have u>i = 2Qq cos 9) . 

The resonance also affects the variance of the differen- 
tial rotation; however, this influence is much weaker and 
disappears at low latitudes (variations in the meridional 
flow have less influence on f2 close to the equator). The 
resonance is not visible in the variance of the four month 
mean flows (panels c, d), which shows the monotonic in- 
crease of the a v and <tq with r c at all latitudes as already 
indicated in Figure ^ 

4. INTERPRETATION 



Our main results presented can be understood in 
terms of different response time scales of meridional flow 
and differential rotation with respect to changes in the 
Reynolds stress. The difference in the response time scale 
follows from the different size of the kinetic energy reser- 
voirs associated with meridional flow and differential ro- 
tation. Here we make a detailed analysis of the energy 
fluxes within the differential rotation model in order to 
estimate typical time scales related to the meridional 
flow and differential rotation. In the following deriva- 
tion we use the assumptions q% <C go, Pi <C po, and 
| S = |V — V a d| <C 1. We also omit terms proportional to 
div(goV m ), which are very small for the low Mach num- 
ber meridional flow and vanish for the stationary solution 
finally conside red. We start from equations (2) to (4) of 
iRempell (J2005) to derive two separate energy equations 
for the meridional flow and differential rotation. To this 
end we multiply equation (2) by gov r and equation (3) 
by Qnvg and add both equations. Rearrangement of the 
different terms leads to: 



d 



+ divF l; 



(r sin 9) 2 g v r< 



grad— 



~~ 2 [RrrE rr + 2R r gE r g + RggEgg + R^E^] 

+v r g g— 

1 

with the energy flux: 

F m = g v m + — - Q v m (rsm6) 2 — {1 

V 2 On J 2 



(G) 



R • V r 



(7) 

Here R denotes the Reynolds stress tensor, E the de- 
forma tions tensor. For the definition we refer t o IRempell 
(2005). Multiplying equation (4) of IRempell l|2005ft bv 

Qo(r sin^) 2 ^ leads to: 



divFr 
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Fig. 3. — Variation of the standard variations of meridional flow 
forcing. Panels (a - b) show the total variance of the time series as 
45° (dotted), and 67.5° (dashed). Panels (c - d) show the variance 
variance caused by inertia oscillations is clearly visible in panels (a, 
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and differential rotation with correlation time scale t c of the random 
function of t c at the top of the domain for the latitudes 22.5° (solid), 
of the 4 month mean for the same variables. The enhancement of the 
b), but not in the four month mean shown in panels (c, d). 



an 

dr 



i on 

r~d6 



-vtQorsmt 



on , 

or 



ldfl. 



+ (r sin 6) 2 g v T 
with the energy flux: 



grad- 



.n 2 



QoV m (fir sin 0) — r sin dflR^ 



(8) 



(9) 



Here R$ denotes a vector with components R r 4> and Re^. 
Since for the boundary conditions we use, the energy 
fluxes [equations Q and |J3}] vanish at the boundaries, 
the divergence terms in equations and JHJ) do not con- 
tribute after integration over the entire domain, leading 
to the energy balances: 



dE m 
dt 

dt 



ic 



= -Qc + Q 
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vise 
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(10) 

(11) 



with the terms 
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dVv t Qor sin 6 



dV = 4tt 
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1 dfl 
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i on 

r~b~0 

A, 



d6r 2 



(12) 



(13) 



denotes the integral over the entire volume of the sphere 
from r = r m i n to r = r max . We emphasize that we 
solve our model only for the northern hemisphere but we 
compute from that the energy conversion for the entire 
sphere. E m and Eq denote energy available in the reser- 
voir of meridional flow and differential rotation, respec- 
tively. For the latter we subtracted the core rotation rate 
to obtain a value representative for the differential rota- 
tion rather than total rotation. Qc denotes the amount 
of energy that is converted by means of the Coriolis force 
(which is a sink for differential rotation and source for 
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Fig. 4. — Schematic view of energy fluxes in differential rotation 
model. 



meridional flow). The terms and represent the 
losses through viscous dissipation for meridional flow and 
differential rotation, respectively. is the energy which 
is converted through the A-effect from internal energy to 
energy of differential rotation, and Q 1 ^ denotes the work 
of the meridional flow against the buoyancy force aris- 
ing from the entropy perturbation within the convection 
zone. 

For the stationary reference model we can write the 
energy balance as: 



Qc 
Qa- 



m 

vise 

n 

vise 



Qc, 



(14) 
(15) 



where the left-hand side terms are the sources and the 
right-hand terms the sinks for the reservoirs of merid- 
ional flow and differential rotation, respectively. Figure^ 
shows a qualitative diagram indicating the energy fluxes 
between the different ene rgy rese r voirs. A quantitative 
analysis of the model 1 of Remnell 1)20051) . which we used 
in this paper as reference model, leads to 



£fo = 1.55 x 10 33 J 
Qa = 6.6 x 10 24 W 



0.017L 



! 



En 
Qa 



Tn = 7^T = 7.5 years , 



= -0.54 Q 
Qc = 0.46 Q\ 



A ■ 



(16) 



© 



for the differential rotation and to 
E m = 10 29 J , 

Qc = 3 x 10 24 J = 0.008 L 

E m 

T m = -pr- = 9.4 hours , 
Qc 

Q™ = -0.02 Q c , 
Q'i? = -0.98 Q c 

for the meridional flow. 

In order to maintain the differential rotation in this 
model, the A -effect has to convert an amount of energy 
equivalent to around 1.7% of the solar luminosity. Com- 
paring this energy flux to the energy stored in the reser- 
voir of differential rotation leads to a time scale of around 



(17) 



8 years for the replenishment of energy. Around 54% of 
the energy converted by the A-effect returns directly to 
the reservoir of internal energy through viscous dissipa- 
tion, while 46% flows through the action of the Coriolis 
force into the reservoir of meridional flow. Dividing the 
energy in this reservoir by the energy flux related to the 
Coriolis force leads to a time scale for the replenishment 
of energy in the meridional flow of only 10 hours, which 
is around a factor of 6500 shorter than the corresponding 
time scale for the differential rotation. The direct viscous 
loss does not play an important role for the meridional 
flow (only around 2% of the energy is dissipated); most 
of the energy returns to the reservoir of internal energy 
through work against the buoyancy force. 

The numbers presented here vary with different model 
assumptions made but the fact that tq is around 4 or- 
ders of magnitude longer than r m is a very robust re- 
sult, since it follows mainly from the different sizes of the 
energy reservoirs of differential rotation and meridional 
flow. Another very robust result is that the viscous dissi- 
pation does not play an important role for the meridional 
flow (unless the assumed turbulent viscosity would be 
more than one order of magnitude larger than the value 
of 5 x 10 8 m 2 s _1 used here). This emphasizes the im- 
portance of the pole-equator entropy variation f or avoid- 
ing th e Taylor-Proudman state, as discussed in iRemnell 
(2005), since the work of the meridional flow against the 
buoyant force is the primary sink of energy. 

From the very short response time scale of r m ~ 10 
hours it can be expected that the meridional flow is 
affected nearly instantaneously by fluctuations in the 
Reynolds stress, while the differential rotation only re- 
sponses to a long term average over several years. As 
a consequence, the relative amplitudes of fluctuations in 
meridional flow and differential rotation are expected to 
scale with a factor of ~ a/to / r m ~ 80, which is seen in 
the results presented before. 

An order of magnitude estimate based on the energy 
yields: 



AE m . 



'V m Av m 

rms rms ' 



Avi 



(18) 
(19) 



Given the fact that the amount of energy transfered into 
the meridional flow is roughly half the energy transferred 
into the differential rotation, the expected change in the 
flow amplitudes scales like 



Av m 1 v n 



rms 



2 v r > 



(20) 



(with v^ ms ~ lkms -1 and v^ ns ~ 5ms" 1 . This is close 
to the ratio we obtained in the more detailed analysis of 
the energy fluxes within the model. 

Since the large kinetic energy reservoir of the differen- 
tial rotation tends to average fluctuations, a systematic 
variation of ctq with correlation time scale and length 
scale is seen in the data (longer lasting, large scale per- 
turbations manifest easier in the differential rotation). 
On the other hand, the meridional flow shows a nearly 
immediate response to changes in the Reynolds stress. 
Therefore a v is mainly determined by the amplitude of 
the random fluctuations in the Reynolds stress assumed 
and not by the correlation time and length scale (unless 
the time scale for fluctuations in the Reynolds stress is 
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much shorter than r m , which is not reasonable in the 
bulk of the convection zone). The only exception is the 
resonance found in case 2, which enhances the amplitude 
by about 50%. 

5. CONCLUSIONS 

We presented a mean field model for differential rota- 
tion and meridional flow including stochastic fluctuations 
in the Reynolds stress driving the differential rotation. 
We found that the variations in the differential rotation 
are around two orders of magnitude less than the vari- 
ability observed in the meridional flow pattern. For ex- 
ample the case 3 discussed above shows a la variation 
of the 4 month mean of less than 0.3% of the rotation 
rate for the differential rotation (in latitudes lower than 
50°), while the la variation for the meridional flow is 
around 8% of the mean flow value (4 month mean) or 
even 20% of the mean flow value if also the short term 
variation is considered. The maximum amplitude of the 
short term meridional flow variability in this models is 
around 10ms -1 and therefore close to 100% of the mean 
flow. 

This result is a consequence of different response time 
scales of the differential rotation and meridional flow to 
changes in the Reynolds stress. Our model shows that 
the energy flows through the reservoir of differential ro- 
tation and meridional flow are comparable, while the ki- 
netic energy of differential rotation exceeds that of merid- 
ional flow by about 4 orders of magnitude. As a conse- 
quence the meridional flow shows a nearly instantaneous 
response (on a time scale of less than a day) to changes 
in the Reynolds stress, while the differential rotation is 
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affected only by changes that have a long term average 
of at least 10 years. 

Our result is in agreement with observations, which in- 
dicate that the la variation of the 4 month differential 
rotation inversions is a few nHz, but also show that the 
variations in the meridional flow can be of up to 10ms -1 . 
However, most of the meridional flow measurements so 
far only address long term variability (time scale of years) 
rather than short term variability, which is the main fo- 
cus of the model presented here. 

We conclude that a fairly significant amount of ran- 
dom noise in the Reynolds stresses driving differential 
rotation can be tolerated without leading to variations 
in the differential rotation contradicting hclioseismic in- 
versions. However, the meridional flow will show signif- 
icant variability on the same order of magnitude as the 
fluctuations of the Reynolds stress. 

We found that the random forcing through the 
Reynolds stress excites inertia oscillations for times scales 
t c ~ 1 - 2% 1 - 4 - 8 days. This resonance enhances 
the total variance a v by about 50% and has a visible but 
much weaker effect on ctq. However, this effect is not vis- 
ible if the variance of the four month mean is considered. 
Therefore this effect could be observable in the variabil- 
ity of daily flow maps of the surface meridional flow but 
not in helioscismic inversions normally considering longer 
averages. 
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